The Structure of TGB C Phases 
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We study the transition from the cholesteric phase to two TGBc phases near the upper critical 
twist k C 2'- the Renn-Lubensky TGBc phase, with layer normal rotating in a plane perpendicular to 
the pitch axis, and the Bordeaux TGBc phase, with the layer normal rotating on a cone parallel to 
the pitch axis. We calculate properties, including order-parameter profiles, of both phases. 



Smectic liquid crystalline phases pj are layered struc- 
tures: they are fluid-like in two-dimensions and solid-like 
in the third. Twist-grain-boundary or TGB phases^, El 
0] are phases of smectic liquid crystals induced by molec- 
ular chirality. They consist of periodically spaced grain 
boundaries, each composed of an array of parallel dis- 
locations, separating smectic slabs as depicted in Figs, 
n and |3 The layer normals N of the slabs rotate in 
discrete jumps across the grain boundaries. These re- 
markable phases are the analog in liquid crystals^] of 
the Abirkosov flux lattice in superconductors [6| with the 
complex smectic mass-density- wave amplitude ip the ana- 
log of the superconducting order parameter, dislocations 
in the grain boundary the analog of vortices, and the 
chiral coupling constant h induced by molecular chirality 
the analog of the external magnetic field H. 

The simplest TGB phase is the TGB^ phase in which 
the smectic slabs between grain boundaries have the char- 
acter of a bulk smectic- A phase in which the layer normal 
N and the director n, specifying the direction of average 
molecular alignment, are parallel to each other in a plane 
perpendicular to the pitch axis along p. In TGBc phases, 
the smectic slabs have the character of a bulk smectic- 
C with n tilted relative to N. Two distinct structures 
for the TGBc phase immediately come to mind. In the 
first 0U, both N and n rotate in the plane perpendicu- 
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FIG. 1: The RL TGBc phase. There is a fixed angle between 
the layer normal N and the director n, which rotate in the 
plane perpendicular to the pitch axis p. 



FIG. 2: The Bordeaux TGBc phase. There is a fixed angle 
between N and n, but n rotates in the plane perpendicular 
to p, and N rotates on a cone whose axis is parallel to p 



lar to p but with a finite angle between them as shown 
in Fig. n We will refer to this as the Renn-Lubensky 
or RL TGBc phase. In the second TGBc phase, first 
discovered 9] and subsequently studied in detail by the 
Bordeaux group and collaborators [Toj . n rotates in the 
plane perpendicular to p, but N rotates on a cone with 
a component parallel to p so that p does not lie parallel 
to the smectic layers as shown in Fig. |5J We will refer to 
this as the Bordeaux or B TGBc phase. No pure form of 
the RL phase has been reported, though phases with two- 
dimensional modulation of the local RL TGBc structure 
have been observed|llj. Though the RL TGBc structure 
may be unstable with respect to these modulations, we 
assume here that it can be stable. We will not discuss 
the TGBc* phaseQ with smectic-C* slabs in which the 
director rotates in a cone from layer to layer. 

Our goal is to study the structure of both the RL and B 
phases near the upper critical field h c i where the TGBc 
phase becomes unstable with respect to the cholesteric 
phase in which the smectic order parameter is zero and 
the director twists in a helical fashion about p with a 
pitch P. Wc follow closely the procedure developed by 
Abrikosov0 in his analysis of the superconducting flux 
phase near the upper critical field H C 2 and applied suc- 
cessfully to the TGBa phase near h C 2 2]. Our analysis 
of the transition to the bordeaux TGB phase is essen- 
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(a) (b) 

FIG. 3: The square amplitude \tp{x, 0, 0)| 2 smectic order pa- 
rameter as a function of x in (a) the Bordeaux and (b) the Rb 
TGBc phases. These figures show the squared slab wavefunc- 
tions \4>{x — nZ;,)| 2 with n — —1 centered at —lb (short dash, 
short space) with n = 1 centered at fo (long dash, short space) 
and with n — (short and long dashes) and \ip(x, 0, 0)| 2 (full 
line). Though |^>(a;)| dies off fairly rapidly, \ip(x, 0, 0)| 2 has a 
robust value at grain-boundary positions x/h = ±0.5. Note 
the asymmetry about x = in the Rb case. 
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TGB A 

Bordeaux 

RL 


-d 2 /du 2 +u 2 
-d 2 /du 2 + u 4 
d 4 /du 4 + u 2 


pl/2^1/2 
p2/3 d l/3 
pl/3 d 2/3 


pl/2^1/2 

P 1/3 d 2/3 

p2/3 d l/3 


pl/2 d l/2 
p2/3 d l/3 
p2/3 d l/3 



TABLE I: 7i u and proportionality of lb, Id, and £ to powers 
of P and d in TGB phases. 



tially identical to that presented by Luk'yanchuk |12| . 
Using a more general model than his, which does not 
permit a stable RL phase, we study both the B and RL 
phases, including their order-parameter profiles, but not 
the TGB2 9 phase he introduced. 

Several results of our analysis are worthy of note. 
The linear stability operator associated with the TGB^ 
phase, like that associated with the Abrikosov phase, 
is a quantum harmonic oscillator Hamiltonian, TL U — 
—d 2 /du 2 + u 2 , where u is a rescaled coordinate along 
p. The same operator associated with the B phase is 
a w 4 -anharmonic oscillator Hamiltonian, ~d 2 jdu 2 + u A , 
whereas that associated with the RL phase is the dual 
of the Bordeaux operator, d 4 /du A + u 2 . The dependence 
of the grain-boundary spacing If,, the dislocation spacing 
Id within a grain boundary, and the smectic coherence 
length £ on the cholesteric pitch P and the layer spacing 
d is different in the three phases as reviewed in Table 
[I] The near equality of If, and Id in the TGByi phase 
and their P 1 / 2 d 1 / 2 dependence on pitch and layer spac- 
ing has been verified 3] . Experimentally ;10], If, is sub- 
stantially larger than Id in the Bordeaux TGBc phase in 
agreement with Tabled Finally, we find that the smectic 
order parameter, though depressed at the grain bound- 
aries, is reasonably constant in the two TGBc phases as 
shown in Fig. [3J We find no evidence in the B phase of 
melted grain boundaries along which ip ~ as suggested 
by Dozov[3. 

To describe the smectic properties of the TGB phases 
near h C 2, we use a slight modification of the Chen- 
Lubensky (CL) model[l4[ for the nematic-to-smectic-A 



and nematic-to-smectic-C transitions. In this model, 
the smectic molecular number density is expressed as 
p — ip + ip* , where ip is the complex mass-density- wave 
amplitude with wavenumbers with magnitude peaked 
near q n = 2n/d. The free energy, — Fh + Fnl is 
the sum of a nonlinear part, Fnl = n9 J d 3 x\ip\ 4 , and a 
part harmonic in ip, 

Fh = y^x[f|Vf + J D|||[Vf|(x)+^]Vf 

+^|[Vi(x)+g \]^| 2 (1) 
+£> NX [Vf|(x) + <$V*[Vi(x) + 9ox]^ + cc], 

where V?,(x) = (n(x) • V) 2 and Vi(x) = V 2 - V 2 |(x). 

With the identification = —C±/2D± and f = r — 
C 2 _/(4:D±), this model is simply a alternative representa- 
tion of the CL model 14] except for the term proportional 
to Z?|| | j_ , which, as we shall see, is needed to stabilize the 
RL phase. When q 2 ± < 0, this model has a phase transi- 
tion from the nematic to the smectic-A phase with ip = 
ipAe lq ° z , where ipA = {-r/g) 1 ^ 2 ■ When qf )± > 0, it has 
transition to the SmC phase with ip = ipce t<ic ' x , where 
ipc = (-r/g) 1 / 2 , with q c = (q ± COS7, q 0± sinj, q ) for 
any angle 7. 

To complete the description of our system, we add the 
Frank free energy, F n , for the director including the con- 
tribution from molecular chirality: 

F n = iydM^i(V-n) 2 + if 2 [n-(Vxn)] 2 

+K 3 [n x (V x n)] 2 + /m- (V x n)}. (2) 

When -0 = 0, the equilibrium state is the cholesteric 
phase determined by F n alone with director 

no(x) = (0, — sin/cox, coskox), (3) 

where k = h/K 2 = 2n/P. 

There are several dimensionless quantities in F = 
Fil> + F n that play a role in our analysis. One is the ratio 
ko/Qo = d/P, which is of order 10~ 2 or less. Our pri- 
mary concern will be the limit in which feo/^o approaches 
zero, and we will consider only leading terms in this ratio. 
Other parameters are the ratios, 

V± = F>±/D\ l: 7711 |_l = D|| ± /L»||, u} = q 2 L0 /q' 2 

and the twist Ginzburg parameter, k 2 = 
(gK 2 /2) 1 / 2 /(4:D\ l q%). In at least one material^. 
771 1 j_ w and r/± <C 1, but there is no a priori reason 
why either of these conditions should not be violated. 
lu = tan a is a measure of the equilibrium tilt angle a 
between n and N. It is more convenient to use the twist 
ko = h/K 2 rather than A as a measure of chirality. The 
critical twist at which the cholesteric phase becomes 
unstable to the TGB phases is k c2 = h c2 /K 2 . 

To determine when the cholesteric phase first becomes 
unstable with respect to the development of smectic or- 
der and to find our variational wave functions 2] for the 
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TGB phases, we calculate the lowest eigenvalues and as- 
sociated eigenf unctions of the harmonic kernel obtained 
from Fh with n(x) replaced by no(x). This kernel K 
is a periodic function of x with period P/2. Its eigen- 
functions are, therefore, plane waves in the yz plane 
that can be expressed as V( x ) = <I> qp (x)e lqpx where 
qp = (0, q y , q z ) and where, as indicated, the form of the 
function $ qp (a;) can depend on qp. When ko = 0, the 
eigenfunctions associated with the lowest eigenvalue of 
K are V'( x ) — e jqc ' x in which qc can have a nonvan- 
ishing x-component. We allow explicitly for this com- 
ponent of ip that varies periodically with x by setting 
^q p (x) = 4>q P (x)e tq!cX and ^(x) = qp (x)e lqx where 
Q = (<fc,qp)- With this form for -0, Fh becomes 



F° H =A 



dx(j>* {x)H(x, 7T, q)(j)q P (x) 



(4) 



where tt — i 1 d/dx is the momentum operator and 



H(x,7r,q) = r + D\\Q\(x,q) + D±Q 2 L (x,TT,q) (5) 
+-D||x[Q|| (x, q)Q±{x,TT, q) + Q±(x,tt, q)Q|| (x, q)], 

where Q||(x,q) = qf { (x) - q% , with gjj(x) = (q • n (x)) 2 , 
and Qj_(x, 7T, q) = TT 2 +2q x % +q 2 ± (x) — g 2 ± , where q\(x) 
T - q 2 (x). 

H(x, 7r, q) is a periodic function of x with a band spec- 
trum and Bloch eigenfunctions. To lowest order in ko/qo, 
however, eigenfunctions are localized at spatial minima 
in r(x, q) = 7i(a;,0, q), which can be approximated by 
the lowest-order terms in a Taylor expansion about these 
minima. For any given q, f (x, q) will have a minimum at 
some x = x m (q). Since fix, q) depends on x only in the 
combination n (a;)-q, x m (U(9)q) = x m (q) + 9/k n , where 
U(9) is the operator that rotates q m through an angle 6 
about p, it is always possible to find a q = q m such that 
i m (q m ) = 0. If 4>(x) is an eigenfunction of TC(x,tt, q m ) 
with energy e, then = 4>(x — #/fco)e lqm( ^' x , where 
q m (fl) = U(0)q m , is an eigenfunction of the harmonic 
kernel of Fh with energy e for all 9. 

Our approach, therefore, is to find those q's that mini- 
mize f (0, q) or, equivalently, those q's for which Qm (0) = 
andQ ± (0,0) = 0. Since q 2 (0) = q 2 z , and q\ (0) = q 2 x + q 2 y , 
it follows that f(0,q) is at its minimum equal to f for 
q z = q and (q x ,q y ) = qo±(cosj, sin 7) for any 7. The 
Bordeaux phase corresponds to 7 = and the RL phase 
to 7 = 7r/2. Having found q m , we can expand Qm and 
Q_l in powers of x and 7r: 



Q|| = -2fc g '7o_L sin7a; + (<? 2 - q 2 ± sin 2 7) + 
Q_L = 7i 2 + 2g i cos77r - Q|| 



(6) 



These expressions simplify in the Bordeaux and RL cases 
to 
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qlk 2 x 2 , Ql = 2q Q± TT - q^x 



2;„2 2 



(7) 



-2fc go9o_La:, 



Q 



PL 



7r 2 - 2fc 9o9oJ_a; (8) 



plus terms, which we show shortly, that yield corrections 
to the lowest order terms in (ko/qo)- If is clear from these 
expressions that the B and RL phases enjoy a sort of dual- 
ity obtained by interchanging x and tt. The Hamiltonian 
for the B phase will have terms proportional to x 4 , tt 2 and 
ttx 2 + x 2 fr, whereas that for the RL phase will have terms 
proportional to tt 4 , x 2 , and xtt 2 + tt 2 x. In the B case, the 
ttx 2 +x 2 fr term can be removed by transforming the wave 
function via 4>b(x) = exp(i/ipx 3 )0p(cc) for an appropri- 



ate choice of /ip while the x t\ 



r 2 term in the RL 



case can be removed by transforming the Fourier trans- 
form (j) RL {k) = j dxer lkx 4> RL (x) to exp(ifj, RL k 3 )4> RL (k). 
In both cases, the eigenfunction 4>b(x) and 4> R l(x) are 
localized near x = over some characteristic length I, 
and it is convenient to express them as functions of the 
unitless variable u — x/l. This leads to the Hamiltonians 
for the RL and B phases expressed to lowest order in u 
and tt u — i~ 1 d/du: 



-i/O 

rLr, 



RL 



q_,0 

n R 



4D||g 4 s 2 w(/co0 2 
D\\qoSi(kol) 4 



1 



TT. 



(fco0 6 4<Z 4 4^" U 
1 4fc 2 77^w 



(k iy 



(9) 



We 



where si = 1 — (f? 2 ^/?7jj and s 2 = 1 + — "^V\\- 
can choose I to make the respective coefficients of tt 4 . and 
7r 2 in H RL and Ti B be unity: 



(kol RL ) 6 = (fc*)/«b)V*l/(4«2W), 
(Mb) 6 - 4(k /q ) 2 {r ]± Lu/ Sl ). 



(10) 



With these choices, H° RL = f + E (T]s 2 ) 1/3 (k /qo) 4/3 [u 2 + 
tt 4 ] and = f + 7 1 2 / 3 Eo(ko/q ) 4/3 [u 4 + TT 2 u }, where £ = 
D||g 4 [16w 2 si] 1 / 3 , are duals to each other with u 2 + tt 4 
and u 4 + tt 2 having the same lowest eigenvalue eo. 

The eigenvalues of both TL RL — r and TC B — f scale as 
Dq 4 (ko/qo) 4 / 3 - Higher order terns in kox and tt neglected 
in Eq. © yield corrections to the dominant (ko/qo) 4 ^ 3 
behavior of both H B — f and H RL — f of order (ko/qo) 2 
or higher. In addition Tt° — f scales as ko/qo when 7 ^ 
0,7r/2, and corrections to H B — f and Ti. RL — f scale, 

respectively, as j 2 (k /qo) 2 ^ 3 and (tt/2 - "f) 2 (k /qo) 2 ^ 3 ■ 
Thus the B and RL phases always have lower energy than 
phases with intermediate values of 7. 

The cholesteric phase becomes unstable at fen — k C 2 
when the smallest eigenvalue of 7i° becomes zero. Thus 
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and near |f| = 0, both k c i(RL) and k C 2(B) scale as \ f | 3 ' 4 . 
Their ratio is k^ 2 /k^ L = (s 2 t 1i _) 1 I 4 . Thus, kf 2 > k* 2 L and 
the cholestreric phase becomes unstable to the B phase 
before the RL phase when s 2 = 1 + — 2??||^ > t]±, 
i.e., when n\± < 1/2, and it becomes unstable to the RL 
phase before the B phase when nn > 1/2. This means 
that the RL phase is not stable in the original CL model 
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FIG. 4: Amplitude of eigenfunctions for (a) the Bordeaux 
and (b) the RL case. Note that the RL wavefunction has an 
oscillatory component and it is slightly asymmetric. 



in which -Dim = 0. Stability of the nematic phase in the 
absence of chirality requires D^D±~D 2 j_ = D^D±s\ > 
or r)?ij_ < r]±. It is clearly possible to satisfy both this 
condition and rju± > 1/2, so that there is a range of 
parameters for which the RL phase is stable. 

The B eigenfunctions are of the form <j)(u) = 
e lfJ -B u ' cf> B (u) where (u 4 + B {u) = e ^ B (u). 
The RL eigenfunctions are of the form 4>rl{u) = 
J(dk/2ir)e~ VRLk 0s(fc), where 4>b(u) and 0s(fc) are 
identical functions of different arguments. 4>b{u) can be 
obtained numerically using the shooting method, and it, 
along with 4>rl(u), is plotted in Fig. 0] The resulting 
eigenvalue is eo = 1.060357.... An excellent approxi- 
mation to 4>b{u) over the entire range of u is 4>(u) = 

exp[-±Au 2 ^Jl + |(w/A) 2 ], where A = 1.035. This func- 
tion satisfies the requirement that 4>{u) — > exp(— -||m 3 |) 
as |u| — * 00. 

To determine the structure of the TGBc phases, we 
construct variational smectic order parameters from the 
degenerate set of lowest eigenfunctions of H(x, w, q): 



where q m (0) = U(6)q m and, as discussed in the para- 
graph preceding Eq. ©, q m has different forms in the 
Bordeaux and RL phases. Following the treatment of the 
cholesteric-to-TGB^ transition[2j, we can write the total 
free energy after minimizing over director fluctuations as 



F 



K 2 q^V 



A 



V2 



(12) 



i>(x) = C^(t>((x - nl b )/l)e 



j?m(e)-x 



(11) 



where (3(lb/l) depends on the separation If, between grain 
boundaries, and A depends on fco an d I but not on 
lb- Thus, the equilibrium value of 4 is that which 
minimizes f3(lb/l), which can be expressed as (['0 4 ]av — 
K 2 _2 ( fc o/go) 2/3 [/]av)/[V' 2 ]av I where / is a complicated 
function of order tp and [g] av — V~ x J cPxg is the spa- 
tial average of g. Carrying out this minimization proce- 
dure using the analytic approximation for 4>b{u), we find 
lb/1 ~ 2.15 for the Bordeaux and lb/1 ~ 3.0 for the RL 
TGBc phases, respectively. To find the spacing between 
dislocations in a grain boundary, we use the geometric 
relation fco = dj (sin a:oWd), where 7r/2 — «o is the angle 
between N and p. Our results for l\, ~ I, ~ Pd/l and 
£ - (Dql/f) 1/2 - d 1/3 P 2/3 are summarized in table [Q 
The wave function tp for our calculated values of lb for 
both TGBc phases are shown in Fig. [3] 

We have presented an overview of the properties of 
the the Bordeaux and RL TGBc phases and the tran- 
sition to them from the cholesteric phase obtained from 
an Abrikosov-like analysis near the upper critical twist 
k C 2 at which the cholesteric phase becomes unstable. In 
a future publication 16], we will present more details of 
our calculations, including a discussion of the transition 
from type I to type II behavior. We will also discuss the 
relation between our work and that of Dozov[T3j. 

We thank Randall Kamien for a careful reading of this 
manuscript. This work was supported in part by the 
National Science Foundation under grant DMR 00-96531. 
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